The strong running coupling at r and Z mass scales from lattice QCD 
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This letter reports on the first computation, from data obtained in lattice QCD with u, d, s and 
c quarks in the sea, of the running strong coupling via the ghost-gluon coupling renormalized in 
the MOM Taylor scheme. We provide with estimates of ajjg"(m^) and Qms(wI) (for which the 
inclusion of the dynamical charm quark makes the running much safer) in very good agreement 
with experimental results. 
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INTRODUCTION 

The confrontation of QCD, the theory for the strong 
interactions, with experiments requires a few inputs: one 
mass parameter for each quark species and the only sur- 
viving parameter in the limit of massless quarks, namely 
Aqcd, the energy scale used as the typical boundary con- 
dition for the integration of the Renormalization Group 
equation for the strong coupling constant, as- Thus, the 
value of the renormalized strong coupling at any scale, 
or equivalently Aqcd, has to be fitted to allow the QCD 
phenomenology to account successfully for experiments 
(see ref. pQ for a description of many precision measure- 
ments of as from different proccesses and at different en- 
ergy scales). The QCD running coupling can be also ob- 
tained from lattice computations, where the lattice spac- 
ing replaces Aqcd as a free parameter to be adjusted 
from experimental numbers. This means that one can 
apply QCD implemented via a lattice-regularized action 
to convert the physical observation used for the lattice 
spacing calibration, as for instance a mass or a decay 
constant, into Aqcd- A review of most of the procedures 
recently implemented to determine the strong coupling 
from the lattice can be found in ref. [5] (see also ref. [3J, 
where we quoted many of the different methods proposed 
in the last few years). 

The present "world average" for the strong coupling 
determinations [4 , usually referred to at the scale of the 
Z° mass, as discussed in ref. pQ, is dominated by the lat- 
tice determination included in the average [5j. Provided 
the importance of a precise and proper knowledge of the 
strong coupling for the LHC cross sections studies and 
its exploration of new physics, independent alternative 
lattice determinations are strongly required. The latter 
is specially true when different lattice actions and proce- 
dures are applied, to gain thus the best possible control 
for any source of systematic uncertainty. Furthermore, 



the current lattice results have been obtained by means of 
simulations including only two degenerate up and down 
sea quarks (N/=2) or, as much (in particular that in [5]), 
also including one more "tuned" to the strange quark 
(N /= 2+l). Now, the European Twisted Mass (ETM) 
collaboration have started a wide-ranging program of lat- 
tice QCD calculation with two light degenerate twisted- 
mass flavours [51 [7] and a heavy doublet for the strange 
and charm dynamical quarks (N/=2+l+l) 0[H]. Within 
this ETM program, we have applied the method to study 
the running of the strong coupling, and so evaluate Aqcd , 
grounded on the lattice determination of the ghost-gluon 
coupling in the so-called MOM Taylor renormalization 
scheme [TUl E]. We are publishing the results of this 
study in two papers: a methodological one [3J, where the 
procedure is dicussed and described in detail and some re- 
sults are presented, and this short letter aimed to update 
and emphasize the phenomenologically relevant results. 
In particular, as far as the lattice gauge fields with 2+1+1 
dynamical flavours which we are exploiting provide with 
a very realistic simulation of QCD at the energy scales 
for the t physics, we are presenting here the estimate for 
the coupling at the r-mass scale and directly comparing 
with the one obtained from r decays. It should be noted 
that including the dynamical charm quark makes also 
safer the running up to the Z°-mass scale. 

THE STRONG COUPLING IN TAYLOR SCHEME 

The starting point for the analysis of this letter shall be 
the strong running coupling renormalized in the MOM- 
like Taylor scheme, 

a T ( M 2 ) = 9 ^fl = lim ^G( M 2 ,A 2 )FV,A 2 ) ,(1) 

47T A->co 47T 

obtained from lattice QCD simulations; where F and G 
stand for the form factors of the two-point ghost and 
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gluon Green functions (dressing functions). The proce- 
dure to compute the coupling defined by ([l}, and from 
it to extract an estimate of A^g, is described in very 
detail in refs. (T01 [TJJ and we recently applied this in 
ref. [3] to compute A^jg from N/ =2+1+1 gauge config- 
urations for several bare couplings, light twisted masses 
and volumes. The prescriptions applied for the appropri- 
ate elimination of discretization artefacts, as the so-called 
i?(4)-extrapolation procedure to cure the artefacts which 
are due to the breaking of the rotational invariance on the 
lattice U2j, needed to be provided with reliable and ex- 
ploitable results were also carefully shown and explained 
in ref. [3 . After this, we are left with the lattice estimates 
of the Taylor coupling computed over a large range of 
momenta that, as can be seen in fig. [TJ can be described 
above around 4 GeV by the following OPE formula: 
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derived from the OPE description of ghost and gluon 
dressing functions [11] . where 7^ can be taken from 
0i] to give for Nf = 4, 1 — jf/p = 27/100, 
and, applying the same method outlined in the ap- 
pendix of ref. [TT], one can take advantage of the 0(a 4 )- 
computations for the Wilson coefficients in ref. [14], and 
obtains R(a, a ) for q Q = 10 GeV (see Eq.(6) of [3 J. The 
purely perturbative running in Eq. (|| is given up to four- 
loops by integration of the ,3- function [3] , where its coef- 
ficients are taken to be defined in Taylor-scheme [Till [TS] . 
The inclusion of the non-perturbative OPE correction is 
unavoidable to properly describe the running of the cou- 
pling over the accesible momentum range, at least above 
around 4 GeV. This however allows to fit both g 2 (A 2 ) 
and At, the Aqcd parameter in Taylor scheme, through 
the comparison of the prediction given by Eq. ([2| and 
the lattice estimate of Taylor coupling. The best-fit of 
Eq. ([2]) to the lattice data published in ref. [3] can be seen 
in the plot of Fig. [T] (dotted line), while the best-fitted 
parameters should be read in Tab. [TJ In this letter, we 
also update the analysis by including an " ad-hoc" correc- 
tion to account for higher power corrections (suggested 
by the detailed analysis of ref. [3]) that allows to extend 
the fitting window down to p ~ 1.7 GeV. Furthermore, 
in addition to the lattice ensembles of gauge configura- 
tions described in ref. [3J, we study 60 more at (3 = 2.1 
(afii — 0.002) and three new ensembles of 50 configura- 
tions at P = 1.9 and am = 0.003, 0.004, 0.005 to peform a 
chiral extrapolation for the ratios of lattice spacings and 
get: a(2.1,0.002)/a(1.9,0) = 0.685(21). Thus, we have 
finally converted our estimates to physical units with the 



previous ratio and the lattice size in the chiral limit at 
P = 1.9 [Hj: 0.08612(42) fm. . 
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FIG. 1: The strong running coupling in Taylor scheme defined 
by obtained over a large momentum range from lattice 
QCD, as described in ref.[3]. The dotted line stands here 
for the best-fit with Eq. H, while the solid one includes a 
higher-order power correction effectively behaving as ~ 1 /p 6 . 



THE WILSON OPE COEFFICIENT AND THE 
HIGHER-POWER CORRECTIONS 

The OPE prediction for ar given by Eq. |2]) is dom- 
inated by the first correction introduced by the non- 
vanishing dimension-two landau-gauge gluon conden- 
sate [TBrl21j . where the Wilson coefficient is applied at 
the 0( a 4 ) -order. In the previous technical paper [3J, we 
provided with a strong indication that the OPE anal- 
ysis is indeed in order: it was clearly shown that the 
lattice data could be only explained by including non- 
perturbative contributions and that the Wilson coeffi- 
cient for the Landau-gauge gluon condensate was needed 
to describe the behaviour of data above p ~ 4 GeV. 

In Fig. [2] the impact of higher-power corrections is 
analized. The plot shows the departure of the lattice 
data for the Taylor coupling from the prediction given 
by Eq. ^ , plotted in terms of the momentum, with log- 
arithmic scales for both axes. The data seem to indi- 
cate that the next-to-leading non-perturbative correction 
is highly dominated by an 1/p 6 term. This might sug- 
gest that the 1/p 4 OPE contributions are negligible when 
compared with the 1 /p 6 ones or that the product of the 
leading 1/p 4 terms and the involved Wilson coefficients 
leave with an effective 1/p 6 behaviour. Anyhow, this im- 
plies that we can effectively describe the lattice data for 
the Taylor coupling with 

c4(p 2 ) = a T (p 2 ) + ± , (3) 

where d is a free parameter to be fitted, which we do 
not attribute to any particular physical meaning, for all 
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FIG. 2: The departure of lattice data from the leading non- 
perturbative OPE prediction for the running coupling plotted 
in logarithmic scales, in terms of the momentum manifestly 
shows a next-to-leading 1/p 6 behaviour; the dotted red line 
stands for the momentum scale, p ~ 1.7 GeV, below which 
the lattice data do not follow any more the 1/p 6 behaviour. 

momenta below p ~ 1.7 GeV. This can be seen in the plot 
of Fig. [T] where the solid black line shows the best-fit of 
Eq. ([3| to the lattice data, with the best-fit parameters 
that can be found in Tab. |T] This table clearly shows 
that the values for the best-fit parameters obtained with 
a formula including the effective higher power correction 
applied for momenta above p ~ 1.7 GeV are compatible 
with the ones obtained without higher power corrections 
and for momenta above p ~ 4 GeV. 



THE STRONG COUPLING IN MS SCHEME 

To obtain the strong MS running coupling from the 
previously estimated values of is rather immediate, as 
the scale-independent AQCD-parameters in both Taylor 
and MS schemes are related through [TT] 

507 - 40iV> 
=e ~ 792 - 487V/ = . 560832 . (4) 

Then, the MS running coupling can be again obtained by 
the integration of /3-function, with the coefficients now in 
Ms-scheme for Ny = 4. Thus, we can apply the estimates 
of Aj^, that can be found in Tab. JTJ to run the coupling 
down to the scale of r mass, below the bottom quark 
mass threshold, and compare the result with the estimate 
from t decays [I], a(m^) = 0.334(14). This can be seen 
graphically in the plot of Fig. [3j while the results will be 
shown in the following. 

The running of a^jg with the two estimates of Aj^g 
from Tab. [TJ and the 1-a error propagation, will produce 
the two following results at the r-mass scale: a^g(m 2 ) = 
0.337(8) and a m (m%) = 0.342(10). If we combine both 
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TABLE I: The parameters for the best-fit of Eq. |2| obtained 
in ref. [3] to lattice data (first row) and the same with Eq. Q 
(second row). The conversion to MS scheme for Aqcd is done 
by applying Eq. Q. The renormalization point for the gluon 
condensate is fixed at fi — 10 GeV. We quote statistical errors 
obtained by applying the Jacknkife's method. 

and conservatively add the errors in quadrature, we will 
be left with the following estimate, 

ais( m ?) = 0.339(13) , (5) 
in very good agreement with the one from r decays. 




q(GeV) 

FIG. 3: The strong MS coupling running for 4 quark flavours 
and for A^jg = 316 MeV (black) and A^g = 324 MeV (blue) 
below the bottom mass threshold. The dotted lines represent 
the one-cr statistical deviations. The red point stand for the 
value of ajyjg^m 2 .) obtained from r decays [T]. 

The determination of a^fg at the Z° mass scale implies 
first to run up to the MS running mass for the bottom 
quark, mb, with /3-coefficients and Aj^g- estimated for 4 
quark flavours, apply next the following matching for- 
mula [1], 

a m =5 ^ = a m =4 ^ ( x + E c «° ( a m =4 ^) ) '( 6 ) 

where the coefficients c„o can be found in ref. |22j . and 
then run from the bottom mass up to the Z° mass 
scale. Thus, from our two estimates of Aj^g, we ob- 
tain: a m (m 2 z ) = 0.1198(9) and a m (m%) = 0.1203(11). 
Again, combining these two results and their errors added 
in quadrature, we will be left with 

«Ms( m D = 0.1200(14) , (7) 
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also in good agreement with the current world aver- 
age [4]: oty^(m 2 z ) = 0.1184(7). It is worth to re- 
mark that this world average is strongly weighted by 
the HPQCD lattice (with 2+1 staggered fermions) re- 
sult [5], oi^{m 2 z ) = 0.1183(8), and by the one from non- 
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singlet structure functions from DIS 



'MS 



0.1142(23), more than 2 er's away from most of the other 
involved estimates. Without these two last estimates, 
the average would be 0.1197(12) and, if one just re- 
places the HPQCD result by the recent one form the 
PACS-CS lattice collaboration (with 2+1 Wilson im- 



a m {m z ) = 0.1205(8)(5), the av- 



proved fermions) [ 
erage would be 0.1190(9), both in still better agreement 
with our estimate. Finally, the world average including 
the present result, with u, d, s and c realistic quarks in 
the loops, instead of the previous lattice one only with 
u,d and s, would be: Q^g(m|) — 0.1191(8). 

It should be noted that we applied two different fitting 
strategies, taking different fitting windows and studying 
the impact of higher order OPE corrections, and no sys- 
tematic effect have been observed. Our error analysis is 
based on the Jacknife method when we account for the 
fitted parameters, while the uncertainties on the lattice 
sizes (and the ratios of them) are properly propagated 
into the final estimates. Some systematic effect related 
to the use of the twisted-mass action for the dynamical 
quarks and to the lattice size determination at the chiral 
limit (not included in our budget error) could also appear 
and can be only excluded by the comparison with other 
lattice and experimental estimates. 



CONCLUSIONS 

We have presented the results for a first computation 
of the strong running coupling from lattice QCD simula- 
tions including w,d,s and c dynamical flavours. We ap- 
plied the procedure of determining the ghost-gluon cou- 
pling renormalized in Taylor scheme over a large mo- 
menta window and then compare this with the running 
obtained in perturbative QCD and improved via non- 
perturbative OPE corrections. That procedure have been 
previously shown to work rather well when analysing lat- 
tice simulations with N/=0 and 2 dynamical flavours and 
so happens here for N/=2+l+l: our estimate for the 
strong running coupling at the r-mass scale nicely agrees 
with those from r-decays and, after being properly prop- 
agated down to the .Z^-mass scale, is pretty consistent 
with most of the estimates applied to obtain the cur- 
rent PDG world average, although slightly larger than 
the N/=2+l lattice result also used for this average. 



[1] S. Bethke, A. H. Hoang, S. Kluth, J. Schieck, I. W. Stew- 
art, et al. (2011), * Temporary entry *, 1110.0016. 
Y. Aoki, PoS LAT2009, 012 (2009), 1005.2339. 

B. Blossier, P. Boucaud, M. Brinet, F. De Soto, X. Du, 
et al. (2011), * Temporary entry *, 1110.5829. 
K. Nakamura et al. (Particle Data Group), J. Phys. G37, 
075021 (2010). 

C. Davies et al. (HPQCD Collaboration), Phys.Rev. 
D78, 114507 (2008), 0807.1687. 

R. Frezzotti, P. A. Grassi, S. Sint, and P. Weisz (Alpha), 
JHEP 08, 058 (2001), hep-lat/0101001. 
R. Frezzotti and G. C. Rossi, Nucl. Phys. Proc. Suppl. 
128, 193 (2004), hep-lat/0311008. 

R. Baron, P. Boucaud, J. Carbonell, A. Deuzeman, 
V. Drach, et al. (ETM Collaboration), JHEP 1006, 111 
(2010), 1004.5284. 

R. Baron, et al. (ETM Collaboration), PoS LAT- 
TICE2010, 123 (2010), 1101.0518. 

P. Boucaud, F. De Soto, J. Leroy, A. Le Yaouanc, 
J. Micheli, O. Pene, and J. Rodriguez-Quintero, 
Phys.Rev. D79, 014508 (2009), 0811.2059. 
B. Blossier et al. (ETM), Phys. Rev. D82, 034510 (2010), 
1005.5290. 

D. Becirevic, et al., Phys. Rev. D60, 094509 (1999), hep- 
ph/9903364; F. de Soto and C. Roiesnel, JHEP 0709, 007 
(2007), 0705.3523. 

J. A. Gracey, Phys. Lett. B552, 101 (2003), hep- 
th/0211144. 

K. Chetyrkin and A. Maier, JHEP 1001, 092 (2010), 
0911.0594. 

K. Chetyrkin and A. Retey (2000), hep-ph/0007088. 
P. Boucaud, A. Le Yaouanc, J. Leroy, J. Micheli, O. Pene, 
and J. Rodriguez-Quintero, Phys.Rev. D63, 114003 
(2001), hep-ph/0101302; Phys.Lett. B493, 315 (2000), 
hep-ph/0008043. 

P. Boucaud et al., JHEP 01, 046 (2002), hep-ph/0107278. 
P. Boucaud, J. Leroy, A. Le Yaouanc, A. Lokhov, 
J. Micheli, et al., JHEP 0601, 037 (2006), hep- 
lat/0507005. 

F. V. Gubarev and V. I. Zakharov, Phys. Lett. B501, 28 
(2001), hep-ph/0010096. 

D. Dudal, R. Sobreiro, S. Sorella, and H. Verschelde, 
Phys.Rev. D72, 014016 (2005), hep-th/0502183. 

E. Ruiz Arriola, P. O. Bowman, and W. Broniowski, 
Phys.Rev. D70, 097505 (2004), hep-ph/0408309. 
K. Chetyrkin, J. H. Kuhn, and C. Sturm, Nucl.Phys. 
B744, 121 (2006), hep-ph/0512060. 

J. Blumlein, H. Bottcher, and A. Guffanti, Nucl.Phys. 
B774, 182 (2007), hep-ph/0607200. 
S. Aoki et al. (PACS-CS), JHEP 10, 053 (2009), 
0906.3906. 



[2 
[3 

K 

[5 

[7: 

[8 
[9 

[io; 

[12 

[13 

[14 

[15 
[16 



[17 
[18 



[19 

[20; 

[21 
[22 
[23 
[24 



Acknowledgements 



We thank the IN2P3 (CNRS-Lyon), IDRIS (CNRS- 
Orsay) and apeNEXT (Rome) computing centers. This 



